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Using the data collected with the D0 detector at yfs = 1.96 TeV, for integrated luminosities 
of about 180 pb _1 , we have measured the ratio of inclusive cross sections for pp — > Z + b jet to 
pp — > Z+jet production. The inclusive Z + b-jet reaction is an important background to searches for 
the Higgs boson in associated ZH production at the Fermilab Tevatron collider. Our measurement is 
the first of its kind, and relies on the Z — > e + e~ and Z — > modes. The combined measurement 

of the ratio yields 0.023 ± 0.005 for hadronic jets with transverse momenta pr > 20 GeV/c and 
pseudorapiditics \g\ < 2.5, consistent with next-to-leading order predictions of the standard model. 



PACS numbers: 14.70.Hp, 14.65.Fy 



Inclusive Z+b-jet production is expected to be a major 
background to Higgs production in the pp — > ZH chan- 
nel, with subsequent Higgs-boson decays into bb. The 
parton-level subprocesses expected to contribute to the 
Z + 6-jet final state are bg — > Zb (where g stands for 



a gluon), and qq — ► Zg, with g — ► bb Q. The process 
bg — > Zb, where the initial b is from the sea of the pro- 
ton parton distribution, is predicted to account for ap- 
proximately two thirds of the total inclusive cross section 
a(pp -> Z + b jet) at y/s = 1.96 TeV. The 6-quark den- 
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sity of the proton influences the production rates of single 
top quarks and the final state hb, with h representing a 
supersymmetric Higgs boson. Consequently, the mea- 
surement of Z + b jet production is an important step in 
constraining the 6-quark density of protons. 

In this Letter, we describe a measurement of the ra- 
tio of production cross sections of inclusive Z + b jets to 
Z+jets. The measurement of the ratio benefits from can- 
cellations of many systematic uncertainties, such as the 
6.5% uncertainty in the luminosity, and therefore allows 
a more precise comparison with theory. 

We search for Z bosons in association with hadronic 
jets in about 180 pb -1 of data collected at the D0 experi- 
ment between August 2002 and September 2003. The D0 
detector at the Fermilab Tevatron collider is a general- 
purpose detector comprising a magnetic central-tracking, 
preshower, calorimeter, and muon systems 0- The 
central-tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker, both located 
within a 2 T superconducting solenoidal magnet. The 
design was optimized for tracking and vertexing capabili- 
ties at pseudorapidities < 3, where r\ = — ln(tan(0/2)) 
and 6 is the polar angle with respect to the proton beam 
direction (z). Particle energies are measured in three 
liquid- argon/uranium calorimeters: A central calorime- 
ter (CC) covers \r)\ < 1.1, and two end calorimeters (EC) 
extend coverage to |nl < 4.2, each calorimeter housed in 
a separate cryostat [3j. Central and forward preshower 
detectors are located just outside of the superconducting 
coil (in front of the calorimetry) , and additional scintilla- 
tors between the CC and EC cryostats provide sampling 
of developing showers for 1.1 < |r/| < 1.4. The muon 
detection system is outside the calorimetry and consists 
of a layer of tracking detectors and scintillation trigger 
counters before 1.8 T iron toroid magnets, followed by 
two similar layers after the toroids. The trigger and data 
acquisition systems arc designed to accommodate high 
luminosities. 

The dielectron sample is selected from the data by 
requiring two clusters of energy in the electromagnetic 
(EM) layers at the trigger level. In the offline selection, 
two EM clusters are each required to have transverse mo- 
mentum pt > 15 GeV/c and \rj\ < 2.5. In addition, 
the shower development in the calorimeter and isolation 
from hadronic activity must be consistent with that ex- 
pected of an electron, and at least one of the EM clusters 
is required to have an associated track to maximize the 
possibility of having a Z boson in the event. The elec- 
tron candidates with matching tracks are to required to 
have a ratio of measured energy in the calorimeter to 
momentum measured with the tracking system consis- 
tent with that expected of an electron. The Z candi- 
dates are selected by requiring a dielectron mass (m ee ) 
of 80 GeV/c 2 < m ee < 100 GeV/c 2 . The Z+jet sam- 
ple is then selected by requiring the presence of at least 
one reconstructed hadronic jet with px > 20 GeV/c and 



\rj\ < 2.5. 

Jets are reconstructed from calorimeter clusters using 
a cone algorithm of cone size A1Z = -J (Ar/) 2 + (A<fi) 2 = 
0.5 in pseudorapidity and azimuth(</>). Hadronic jets are 
required to have an associated cluster of tracks ("track 
jets"). This requirement reduces background from noise 
in the calorimeter. Track jets are found by applying 
a cone track clustering algorithm of size A1Z = 0.5 
with a seed track of px > 1-0 GeV/c, to tracks of 
Pt > 0.5 GeV/c that are close to the primary interac- 
tion vertex (whose determination is discussed below). A 
track jet can consist of two or more tracks. 

A "taggable" jet is a calorimeter jet with a matching 
track jet within A1Z < 0.5. Applying the taggability cri- 
terion to 2,661 jets in 2,219 Z(ee) candidate events in the 
mass Z window yields 1,658 events. Based on side bands 
to the Z mass window, 121 ±4 events are estimated to be 
from background sources. The main background is from 
multijet production where two jets mimic EM objects, 
with one of the objects having an overlapping track that 
passes the track-matching criteria. The taggability per 
jet is (75 ± 1)% after background subtraction. 

The dimuon sample is defined by the detection of at 
least one muon candidate at the trigger level. In the off- 
line selection, two isolated muons are required to be of 
opposite charge, and to have px > 15 GeV/c and \-q\ < 2 
with trajectories in the muon spectrometer matched to 
tracks in the central-tracking detector. Muon isolation is 
based on the transverse component of the muon momen- 
tum relative to the combined momenta of muon and the 
closest calorimeter jet in {rj,4>) space, and is pTrei > 10 
GeV/c. The Z candidates are selected by requiring a 
dimuon mass of 65 GeV/c 2 < m w < 115 GeV/c 2 . The 
Z mass window is larger than in the dielectron channel 
due to worse momentum resolution for high px muons. 
The criteria for reconstructed hadronic jets are the same 
as in the dielectron channel. A total of 1,406 events re- 
main after the requirement that there be at least one tag- 
gable jet. The main background in this channel is from bb 
production, where both b jets contain muons that satisfy 
the isolation criterion (referred to as bb background) . The 
isolation efficiencies of muons from Z and bb are expected 
to be different, since, for the latter, a hadronic jet would 
be expected to be close to the muon. By performing 
fits to dimuon mass spectra, where the background con- 
tributes to the continuum, samples with different number 
of isolated muons are analyzed to measure the isolation 
efficiencies and background rates. From such analyses, 
we estimate that the background contribution to the fi- 
nal sample with two isolated muons is 17.5 ± 4.1 events. 

Figure ^ shows distributions in transverse momentum 
of taggable jets for both channels (points with error bars), 
compared to a Z+jet Monte Carlo (MC) generated with 
ALPGEN 0, using PYTHIA || for parton showering and 
hadronization. Also shown is a background estimation 
based on data obtained from samples that are in the 
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FIG. 1: The pr distribution of taggable jets in dielectron and 
dimuon channels compared to Z+jet ALPGEN with pythia 
showering and full detector simulation (open histogram), and 
background (multijet for ee channel and bb background for fi^j, 
channel) derived from data. The error bars on the data points 
are statistical. The prediction is normalized to the data, as 
described in the text. 



side band for the dielectron channel, or fail the isola- 
tion criterion for the dimuon channel. The background 
distribution is normalized to the number of background 
events estimated in the selected sample. The simulated 
signal is then normalized so that the total agrees with 
the measurement in Fig. ^ Within the uncertainty of 
the jet energy scale (JES), indicated by the darker shad- 
ing about the expectation, the shape of the distribution 
is well-described by the simulation. 

The b quark fragments into a B hadron, which is iden- 
tified by a displaced secondary vertex that is separated 
from the primary vertex. The reconstruction of sec- 
ondary vertices proceeds in two steps. First, the pri- 
mary interaction vertex is identified, and then additional 
nearby vertices are reconstructed. In the high luminos- 
ity environment of the Fermilab Tevatron collider, there 
can be more than one interaction per beam crossing, one 
of which is likely to have triggered the recorded event. 
The interaction region in D0 has a root mean-square 
width of w 25 cm along z, with a transverse beam size 
of ~ 30 fim. It is possible to distinguish the main hard- 
interaction vertex from any additional soft interactions 
because the vertices are normally well-separated along 
z. Primary interaction vertices are reconstructed in two 
passes. In the first pass, all tracks present in an event 
are used to find seed vertices using an iterative method, 
where tracks that contribute to a fit to a common vertex 
with a x 2 /d.o.f. greater than some chosen threshold are 
removed. The fit is repeated until a stable set of seeds 
is obtained. The seed vertices are then used in a second 



pass to fit all tracks within a certain distance-of-closest- 
approach to any seed. This improves the position resolu- 
tion on the vertex, since the fit is less affected by poorly 
reconstructed tracks. The px distribution of the associ- 
ated tracks is then used to select the primary interaction 
vertex (PV). 

A 6-jet tagging algorithm for secondary vertices (SV) is 
used to identify heavy-quark jets in the analysis. Tracks 
that are displaced from the PV in the transverse plane 
are used as seeds to find secondary vertices. First, a 
fixed-cone jet algorithm of A1Z = 0.5 is used to cluster 
the tracks to form track-jets. Tracks are required to have 
hits in at least two layers of the SMT, px > 0.5 GeV/c, 
and be within 0.15 cm in the plane transverse to z and 
0.40 cm in z relative to the PV. Tracks identified as aris- 
ing from Kg and A decays or photon conversions are not 
considered. Any pair of tracks within a track-jet with an 
impact parameter relative to the hard-interaction vertex 
(distance of closest approach - dca - of a track to a ver- 
tex in the plane transverse to the z direction) divided by 
their uncertainties (<7d ca ), dca/odca > 3 is used as a seed 
for secondary vertices. Additional tracks are attached it- 
eratively to the seed vertices if their ^-contribution to 
the vertex fit is consistent with originating from the ver- 
tex. A secondary vertex consists of two or more tracks. 
The momentum vector of the SV is defined as the vector 
sum of track momenta. Finally, good-quality secondary 
vertices are selected based on the decay length (distance 
between PV and SV), collinearity of the vertex momen- 
tum with the direction from PV to SV, and vertex-fit 
X 2 - A jet is considered b tagged when it is taggable and 
has at least one secondary vertex, with a decay-length 
transverse to the PV {L xy ) divided by its uncertainty 
Lxy/o'xy > 7, associated with it. A secondary vertex is 
associated to a jet if the opening angle between the direc- 
tion of the calorimeter-based jet axis and the momentum 
vector of the SV is ATI < 0.5. 

The ^-tagging efficiency (e&) and the light-flavor tag- 
ging rate (el) of the fe-tagging algorithm are parametrized 
as functions of jet px and 77. The parametrization of eb is 
derived from a different data sample using events with 
jets containing muons (muonic jets), which are domi- 
nated by b jets, but also have contributions from light 
quark jets, gluon jets, and charm jets. The 6-tagging 
efficiency is extracted from the heavy-flavor component 
in this muonic jet sample. The light-flavor tagging rate 
is also derived from data, after compensating for effects 
of displaced vertices that do not originate from heavy- 
flavor decay (Kg, A, and photon conversions). Different 
types of samples are used to determine ez, and Cb, and 
the spreads are taken as systematic uncertainties. 

A comparison of inclusive Z+jet events, generated 
with the ALPGEN leading-order matrix element and 
pythia for showering, with inclusive Z + b events gener- 
ated with pythia, shows good agreement for jet px and 
r\ distributions. We therefore use the shapes of px and 
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■q derived from the Z+jet data sample to estimate the 
expected 6-tagging efficiency and the light-flavor tagging 
("mistag") rate. The average 6-tagging efficiency and 
mistag rate per jet, averaged over p? and 77, are found to 
be (32.8 ±1. 3)% and (0.25 ± 0.02)%, respectively, for the 
dielectron channel. Corresponding values for the dimuon 
channel are (33. f ± 1.1)% and (0.24 ± 0.02)%. To ob- 
tain the event mistag rate, we take into consideration jet 
multiplicity, and measure the event mistag rate of 0.28% 
(0.27%) for the dielectron (dimuon) channel. 

Since tb is derived from events with a muon embedded 
in a jet, whereas most of the 6-tagged jets do not contain 
such muons, the difference in 6-tagging efficiencies for 
hadronic 6 jets and muonic 6 jets is derived from MC, 
and the ratio is used to correct e b . We cannot at this 
point derive the charm tagging efficiency (e c ) from data, 
so we rely on pythia MC to compare Z — ► 66 and Z — > 
cc samples. We assume that (e c /e b ) data = (e c /e fc ) M c = 
0.266 ±0.003. 

The jet taggability, tz,, is measured using data to be 
(75 ± 1)%, while that for b jets, t b , is obtained from MC, 
and scaled such that {t b )d&t a = (*i)data X (|^)mc. The 
result is (tbjdata — (79.2±1.3)% for the dielectron channel 
and (80.7 ± 1.1)% for the dimuon channel. We assume 
that the taggability of charm jets is same as tb- 

After applying b tagging, 27 Z(— > ee) + 6-jet candidate 
events are left, with an expected background from the 
Drell-Yan ee continuum and multijet background of 4.2 ± 
1.4 events based on the side-band subtraction method. In 
the dimuon channel, 22 events are observed with 5.0 ± 1.1 
events from bb background. 

After subtracting the background contributions, two 
equations, one before and the other after the requirement 
of b tagging, determine the contributions from different 
flavors in the remaining events: 

AWoreb-tag = t b N b + t c N c + t L N L (1) 

iVb-taggcd = e b t b N b + e c t c N c + e L t L N L , (2) 

where Nb, N c and Nl are the numbers of events with 
6, c and light jets, respectively; U are the taggabilitics 
per event for different jet types; and the q are the cor- 
responding mean event-tagging efficiencies. We assume 
that the tagging efficiencies per jet, lb and e c , are the 
same as the tagging efficiencies per event. Equations 
and J5J) have three unknowns. We take the theoretical 
prediction of N c = 1.69iV& l], to provide a solution of 
Eqs. and @ for N bl N c and N L . 

The ratio cr(pp — ► Z + b jet) / 'a(pp — ► Z+jct) = 
N b /(N b + N C + N L ) is found to be 0.026 ± 0.007 for the 
dielectron channel and 0.020±0.005 for the dimuon chan- 
nel, where the errors are purely statistical. The combined 
ratio, using the statistical weighting of the number of ob- 
served Z+]d candidates, is 0.023 ± 0.004. The shape of 
the pt spectrum for 6-tagged jets and the significance of 
decay lengths of secondary vertices are compared to the 
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FIG. 2: (a) The pr spectrum for 6-tagged jets, (b) Distribu- 
tion in decay-length significance of secondary vertices in the 
transverse plane, without the requirement on decay- length 
significance. All error bars are statistical. 



TABLE I: Systematic uncertainties for the combined ratio of 
cross sections, showing the impact of ±1 standard deviation 
changes in contributions. 



Source 


Upward (%) Downward (%) 


Jet energy scale 


5.7 


6.7 


Background estimate 


5.6 


5.3 


Z + (QQ) 


0.0 


5.5 


Mistag rate 


3.5 


3.3 


6/c tagging efficiency 


3.1 


3.2 


Taggability 


1.9 


3.7 


Correction for hadronic jet 


2.0 


1.7 


Jet reconstruction efficiency 


1.8 


1.8 


a(Z + c)/a{Z + b) 


2.8 


2.9 


Total (added in quadrature) 


10.2 


12.3 



sum of background and Z + b MC in Fig. |3 The con- 
tribution of each component is given by the solution to 
Eqs. and J2J. The distribution of the decay- length 
significance for secondary vertices shows clear evidence 
for a heavy-flavor component in the 6-tagged candidate 
events. 

Sources of systematic uncertainty in the ratio include: 

i) Jet energy scale. The JES is varied within its un- 
certainty. The JES for hadronic 6 jets is assumed to 
be the same as that for light-flavored jets, whereas 
in MC some differences are observed, and this effect 
is included as part of the JES uncertainty. 

ii) Different methods of estimating background. The 
background is varied by its measured uncertainty 
and the ratio is recalculated. 

iii) Jets that contain a 66 or cc pair from gluon split- 
ting. These jets have a higher tagging probability. 
The expected contribution is taken from theory [jj] , 
and the relative increase in 6(c)-tagging efficiency 
is estimated from MC. This is labeled as Z+ (QQ) 
in Table H| 
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iv) Mistag rate for light jets, which depends on the 
type of jet sample. Using events collected from 
hadronic jet triggers, the light-jet tagging efficiency 
is measured to be 0.23%, and for a sample of events 
with an enhanced EM fraction and small imbalance 
in overall pr, this is 0.26%. A tagging efficiency of 
0.25% per jet (or 0.28% per event) is obtained for 
the combined data. 

v) Uncertainty in tagging efficiency for b and c jets is 
obtained by varying the efficiency by a ±1 stan- 
dard deviation, assuming complete correlation in 
the ratio of extracted cross sections. Also, for c 
jets, there is additional uncertainty from the e c /eb 
ratio obtained from MC. ei is varied, as above, to 
estimate this effect. 

vi) A small difference observed in for different 
MC samples of Z + b jet/Z+light jet, and Z -> 
bb/Z — > qq is taken into account. 

vii) Differences in tagging efficiency between hadronic 
jets and those containing muons. The 6-tagging 
efficiency is measured in data using muonic jets. 
The tagging efficiency for hadronic jets is estimated 
to be 86% of that of muonic jets, as derived from 
Z bb MC. The same ratio in Z + bb MC is mea- 
sured to be 84%, and the difference of 2% is taken 
as a systematic uncertainty. 

viii) Different p^-dependence in jet reconstruction for 
light, b, and c jets, measured using MC samples, is 
accounted for as a systematic uncertainty. 

ix) Uncertainty from theory for the ratio a(Z + 
c jet)/a(Z + b jet) = N c /Nb is estimated as 9.5% 


The effects of systematic uncertainties on the com- 
bined measurement are listed in Table [I] All these uncer- 
tainties are assumed to be completely correlated for the 
two channels, except for that due to background estima- 
tion. Folding these uncertainties together, yields a ratio 
of 0.023 ± 0.004(stat)±°;°^(syst). This measurement is 
in good agreement with the next-to-leading order (NLO) 
prediction of 0.018 ± 0.004 @. 

In summary, we have presented the first inclusive mea- 
surement of 6-jet production in association with Z bosons 



at the Fermilab Tevatron collider, which is a background 
to the standard-model Higgs searches in the ZH pro- 
duction channel. The measurement is in agreement with 
the NLO calculations and can be used to constrain the 
6-quark density of proton. 
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